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We analyze the steady-state transmission of high-momentum (high-/c) electromagnetic waves 
through metal-semiconductor multilayer systems with loss and gain in the near-infrared (NIR). 

Using a semi-classical optical gain model in conjunction with the scattering matrix method (SMM), 
we study indium gallium arsenide phosphide (InGaAsP) quantum wells as the active semiconductor, 
in combination with the metals, aluminum-doped zinc oxide (AZO) and silver (Ag). Under mod¬ 
erate external pumping levels, we find that NIR transmission through Ag/InGaAsP systems may 
be enhanced by several orders of magnitude relative to the unpumped case, over a large angular 
and frequency bandwidth. Conversely, transmission enhancement through AZO/InGaAsP systems 
is orders of magnitude smaller, and has a strong frequency dependence. We discuss the relative 
importance of Purcell enhancement on our results and validate analytical calculations based on the 
SMM with numerical finite-difference time domain simulations. 


I. INTRODUCTION 

Metal-dielectric (MD) interfaces support optical sur¬ 
face waves whose effective wavelength is less than the 
smallest wavelength achievable in dielectrics 1 . The cou¬ 
pling of multiple MD interfaces enables even smaller ef¬ 
fective wavelengths, or conversely higher effective indices, 
or so-called high-fc propagation, of the optical waves*®. 
The phenomenon of large effective indices and propa¬ 
gating high-& modes lies at the heart of sub-diffract ion- 
limited imaging device!^ 7 , sensors based on plasmonic 
resonances^, as well as applications of hyperbolic meta¬ 
materials (HMMs), including asymmetric transmission 
device^, nonlinear opticd^J and lifetime reduction of 
dye molecule^ 12 13 and quantum dots 1 ®. Inherent to the 
increase in the effective index in MD syste ms, is a con¬ 
comitant increase in dissipation losses 15 16 . The tradeoff 
between optical confinement and losses is a hallmark of 
plasmonic and hyperbolic media, and may be considered 
one of the greatest challenges to the w idespr ead realiza¬ 
tion of integrated plasmonic technology 17 18 . 

Strategies for reducing losses in plasmonic and HMMs 
have therefore been proposed, including clever design of 
passive structures and incorporation of active, i.e. gain 

medisP^JHsj 

An exemplar y cas e of the former strategy 
is the fishnet metamaterialwhereby air voids in an 
otherwise continuous MD multilayer open new channels 
for optical transmission. Even in these structures though, 
gain media are still necessary to obtain low-loss trans¬ 
mission at optical frequencies, as was demonstrated via 
infiltrated optical dyes in the core of the structure 25 . 

While dye molecules are straightforward to model 
as non-interacting two-level systems and convenient for 
proof-of-concept experiments^®^, an ideal gain media in 
plasmonic systems would be electronically addressable 29 . 
In this case, semiconductor heterostructures become 
strong candidates for incorporation into plasmonic sys¬ 
tems. III-V semiconductors emitting in the near-infrared 


(NIR) have been studied in conjunction with surface 
plasmon amplificatiorP^, and successfully implemented in 
sub-wavelength metal-coated and plasmonic sources 30 32 . 
However, the incorporation of active III-V materials in 
multilayer MD systems with hyperbolic dispersion re¬ 
mains to be demonstrated in the NIR. 

In this work, we apply the scattering matrix method— 
(SMM), to analyze the possibility for using indium gal¬ 
lium arsenide phosphide (InGaAsP) multiple quantum 
wells (MQWs) as a gain media in NIR MD systems. In¬ 
GaAsP is a mature quaternary III-V compound that is 
epitaxially grown on indium phosphide (InP) substrates. 
This immediately poses challenges for its incorporation 
in a multilayer MD system. We envision deposition of 
metallic thin films into finely patterned trenches between 
InGaAsP pillars formed by nanoimprint lithography and 
subsequent reactive ion etching. We assume that the 
smallest layer thicknesses enabled by this process are 
30nm. The resulting MD system, shown schematically 
in Fig{l] supports volume plasmon polaritond^ (VPPs) 
with large effective indices and large loss in the absence 
of external pumping, the hallmark of strongly coupled 
MD interfaces. These high-& VPPs propagate normal to 
the epitaxial growth direction, which is a novel multilayer 
configuration, potentially more suitable to waveguide- 
integrated HMMs, compared to the more common large- 
area multilayers which are conformal to the wafer sub¬ 
strate. In the following, we show that under moderate 
external pumping conditions, MD systems composed of 
InGaAsP MQWs may support transmission of high-index 
modes enhanced by several orders of magnitude, relative 
to the unpumped case. 

In Section |H| we introduce the optical gain model em¬ 
ployed, in conjunction with the SMM. We present the 
main results of the report in Section |III[ namely, the 
NIR transmission characteristics of multilayer MD sys¬ 
tems based on InGaAsP MQWs. We investigate In¬ 
GaAsP in combination with aluminum-doped zinc ox- 
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FIG. 1. Schematic of multilayer metal-InGaAsP MQW 
system. The heterostructure, with t we ii = 10nm and 
tbarrier— 20nm is grown in the y-direction and TM-polarized 
light propagates in the z-direction. 


ide (AZO), which is considered a low-loss transparent 
conducting oxide and alternative plasmonic material for 
NIR application^®^, and in combination with silver, 
which is well known as th e n oble metal with highest 
conductivity 37 . In Section IV, we discuss the relative 
importance of Purcell enhancement on our results. We 
also validate the analytical SMM results with numerical 
finite-difference time-domain (FDTD) simulations. We 
then discuss the limitations of our results imposed by 
the assumptions of our models. We conclude the report 
in Section lYl 


II. METHODS 


The quaternary III-V compound that we consider is 
an ln x Gdii- x AsyPi-y MQW system where (x=0.564, 
y=0.933) and (x=0.737, y=0.569) for the wells and 
barriers, respectively. The nominal room-temperature 
bandgap energies of the barrier and well materials are 
0.953eV (AG=1.3/im) and 0.774eV (AG=1.6/im), respec¬ 
tively. This system can be epitaxially grown on an InP 
substrate and has been experimentally utili zed i n near- 
infrared subwavelength semiconductor laser®® 

In the dipole-approximation the peak magnitude of the 
optical gain is governed by the transition matrix element, 
Mt . Because MD multilayers only support propaga¬ 
tion of TM-polarized modest, we consider TM-polarized 
light emission from the quantum wells. To justify the 
use of a steady-state analysis and the non-inclusion of 
dynamic nonlinear effects, and for simplicity, we con¬ 
cern ourselves with moderate pumping levels, about two 
times above the transparency inversion density. Conse¬ 
quently, we may assume that most of the transitions oc¬ 
cur at the band-edge such that recombination between 
electrons and holes with vanishing transverse momen¬ 
tum dominate®! In this limit, only transitions from 
the conduction band to the light-hole band contribute 
to the matrix element for TM-polarization 39 ! We ap¬ 
proximate the bandstructure of InGaAsP with parabolic 
conduction and valence bands, characterized by the con¬ 


duction and light-hole effective masses of 0.0481rao and 
0.0537mo, respectivel}®, where rao=9.1xlO -31 lkg is the 
free electron m ass. The magnitude of the matrix element 
used is given by * 41 1 42 1 |Mr| 2 = (24.9eV)rao/3, which corre¬ 
sponds to a dipole length of |x| ^2.8nm at Ao=1500nm. 

Employing this gain model we obtain complex permit¬ 
tivity values for a lOnm InGaAsP QW at carrier densi¬ 
ties representative of absorbing (Abs) and moderately in¬ 
verted (Inv) states, N= 1.0x16cm -3 and N= 5.0x18cm -3 , 
respectively. At Ao=1500nm, the imaginary parts of 
these values are ep^Abs ”=rg(7V=1.0xl6cm 3 )=±0.139 
and eD ? / nv ”=rg(V=5.0xl8cm _3 )=-0.143, where T= 1/3 
is the relative area of wells in the InGaAsP heterostruc¬ 
ture. Using Kramers-Kronig relation®*, the real parts 
at A 0 =1500nm are e' D , Abs =11.997 and e' D Inv =ll.SSl, re¬ 
spectively. For the sake of simplicity, we used a carrier- 
independent but frequency dependent real permittivity, 
which is 6^=11.914 at Ao=1500nm, and justified because 
the transmission properties studied are dominated by loss 
and gain. Relating the imaginary permittivity to the 
linear gain coefficient^, g=-&oFD”(e^) -1 / 2 , the value of 
e£>” = ±0.14 corresponds to a loss/gain per unit length 
of ±1700cm -1 at 1500nm. The transparency (Tra) 
conditior®, eD,Tra”= 0? for InGaAsP at Ao=1500nm cor¬ 
responds to about N= 2.0x18cm -3 . The complex per¬ 
mittivity values for AZO® and A^SS at Ao=1500nm are 
-0.392±i0.139 and -122.190H-i3.115, respectively. Addi¬ 
tional de tails on the gain model may be found in the 

reference d 39 * 45 !. 

Combined with the presented optical gain model, we 
use the SMM as implemented by Krayzel et al® to study 
the transmission properties of an idealization of the struc¬ 
ture of Fig|l] In comparison to simpler methods such as 
effective medium theory 33 , Bloch’s theorem®, and the 
transfer matrix method 20 2 ®Zl, the SMM is well known as 
a superior method for analyzing systems with evanescent 
wave®EU and strong loss/gair®. The idealized struc¬ 
ture includes loss and gain, is infinite in the transverse 
(x-y) plane and surrounded by a uniform medium in the 
longitudinal (z) direction. To couple to the high-& modes 
supported by the MD multilayer of FigJTJ the surround¬ 
ing medium takes the form of a prism. The prism behaves 
as a numerical simplificatio n to the physically equivalent 
grating-coupling technique 10 12 I 14 I 51 I 52 J which is required 
for practical excitation of the high-fc modes supported by 
each system. 

To determine a suitable value of the prism permittivity, 
6p, we first determine the bounds of the high-fc trans¬ 
mission windows in each MD system. This is achieved 
using Bloch’s theorem® in the absence of losses by lo¬ 
cating the regions in momentum space wherein the real 
part of the Bloch vector is purely real. For MD systems 
with 30nm layer thicknesses and Ao=1500nm, Fig§a) 
and[2|b) show the real (solid blue curve) and imaginary 
(dotted blue curve) parts of the Bloch vector, indicating 
that transmission windows exist when 0<k x <2.83ko and 
4.91fco<fc x <6.56fco in AZO/InGaAsP and Ag/InGaAsP 
systems, respectively. In the limit of effective medium 
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30/30nm AZO/InGaAsP, X 0 =1500nm lDJ 30/30nm Ag/InGaAsP, A. 0 =1500nm 




FIG. 2. Wave-vector diagrams for 30/30nm (a) 

AZO/InGaAsP and (b) Ag/InGaAsP systems at 1500nm 
with losses omitted. (solid blue line=real part of Bloch 
vector, dotted blue line=imaginary part of Bloch vector, 
dashed red line=effective medium theory prediction of real 
part of Bloch vector) The transmission windows in (a) 
and (b) extend from 0<k x <2.83ko and 4.91ko<k x <6.56ko , 
respectively. 


theory (dashed red curves), Fig. §a) and [2|b) are rep¬ 
resentative of type-I and type-II HMMs, respectively. 
Based on these high-A: windows, we set ep=64 for the 
SMM calculations. 


III. RESULTS 

Transmission and reflection at Ao=1500nm through 
a 10-period AZO/InGaAsP structure with 30nm layer 
thicknesses, is shown as a function of incident angle 
in Figj3^a) and lb), respectively. The resonances 
extending from normal incidence to 9i nc = 20° trans¬ 
late into a high-/c window that extends to roughly 
k x =8 sin(20°)~2.8&o, consistent with predictions based 
on Bloch’s theorem (Fig. [2ja)). The effect of inverting 
the carrier population in absolute terms appears most 
pronounced near normal incidence and becomes weaker 
as the incident angle increases. The strong absorption in 
this system is apparent in the local reflection minimum 
around 0i nc =7. 5°. The lack of corresponding transmis¬ 
sion maximum indicates strong dissipation, even when 
the semiconductor is inverted. 

Transmission and reflection through the 
AZO/InGaAsP system is best appreciated when 
contrasted to the behavior of the Ag/InGaAsP system. 
In Fig(3jc) and |3jd), transmission and reflection for the 
latter system are shown, with identical conditions to the 
AZO-based structure. The transmission window spans 
approximately 5-Oko<k x <6.5ko, again consistent with 
Bloch’s theorem (Fig. |2jb)). The effect of inverting the 
carrier population becomes dramatic in both absolute 
and relative terms. For several resonances between 5.0/co 
and 5.7fco, the transmission and reflection exceed unity, 
indicating active behavior of the system. 

To quantify the effect of the inverted carrier population 
on the behavior of the MD systems, we define the relative 
transmission enhancement factor ( TEF ) as 

TEF = T(5xl0 18 cm“ 3 )/T(lxl0 16 cm“ 3 ), (1) 


lOp 30/30nm AZO/InGaAsP, s p =64, ?„ 0 =1500nm 


lOp 30/30nm AZO/InGaAsP, s p =64, *. 0 =1500nm 




lOp 30/30nm Ag/InGaAsP, s p =64, X Q =1500nm 


lOp 30/30nm Ag/InGaAsP, s p =64, A 0 =1500nm 




FIG. 3. (a,c) Transmission and (b,d) reflection for TM- 

polarized light of wavelength 1500nm, incident on 10- 
period (a,b) AZO/InGaAsP and (c,d) Ag/InGaAsP multi¬ 
layer with 30nm layers, coupled via prism with ep=64. (Abs: 
iV=lxl0 16 cm -3 ; Tra: ep”=0; Inv: iV=5xl0 18 cm -3 ; [0]: 
cm” =cd” =0) 


the transmission under inversion with respect to trans¬ 
mission under absorption of the InGaAsP MQW. 
FigQa) shows the TEF for the AZO-based multilayer 
at several different incident wavelengths. This system ex¬ 
hibits significant dispersion due to the proximity of the 
plasma frequency of AZO to the NIR. For shorter wave¬ 
lengths, enhancement is negligible, indicating that real¬ 
istic gain levels cannot compete with the strong damping 
at the plasmon resonance. As the wavelength increases 
to 1550nm, a modest TEF appears, increasing with the 
incident angle. However, this relative enhancement must 
be celebrated cautiously. The improvement in the TEF 
with increasing angle of incidence is offset by a drop in 
the absolute transmission (Figj3ja)). Consequently, the 
overall effect of carrier inversion is quite small in the 
AZO/InGaAsP system. 

Again, these results are best appreciated by contrast¬ 
ing to the Ag/InGaAsP system, the TEF of which is 
shown in Fig®b). Relative enhancements close to a fac¬ 
tor of 10 are observed over the range of incident angles 
supported by this system, with several prominent peaks 
exceeding factors of 100 near the edges of the transmis¬ 
sion window. In comparison to the AZO/InGaAsP sys¬ 
tem, the TEF has weak dispersion and is consistently 
larger in magnitude, indicating that the effect of the 
inversion in the Ag/InGaAsP system is much stronger. 
Thus, if efficient coupling to this system is achieved, tun¬ 
able and/or active, pump-dependent behavior may be 
feasible, enabling, for example, extremely compact and 
electronically addressable optical amplifiers, sources, and 
mixers. 














































































4 






FIG. 4. Transmission enhancement factor ( TEF ) for TM- 
polarized light incident on 10-period (a) AZO/InGaAsP and 
(b) Ag/InGaAsP multilayer with 30nm layers, coupled via 
prism with ep=64. 


IV. DISCUSSION 


Prior work on gain-compensated plasmonic systems 
has emp hasize d that the Purcell effect has a delete¬ 
rious rold- 29 47 . Specifically, while the carrier densi¬ 
ties required for significant improvement in transmis¬ 
sion appear practically achievable, the required current 
densities may be extremely high 29 . This is a conse¬ 
quence of the large optical density of states in plasmonic 
and HMMs and the related high rates of spontaneous 
emissiorF^. We stress, however, that the Purcell effect 
in MD systems is most significant near the plasma fre¬ 
quency of the constituent metal and becomes less impor¬ 
tant as the frequency decreases. Using the expression 29 
Fp=l-\-'KTk z ,]jk x uj(dk x / duo) / (ejjko) 2, for the Purcell fac¬ 
tor, Fp, at a single MD interface, where k z ^ is the lon¬ 
gitudinal wave-vector component in the dielectric, and 
uj and ko are the angular frequency and vacuum wave- 
number, respectively, we immediately notice that the 
Purcell enhancement in the AZO/InGaAsP system is sig¬ 
nificantly greater than unity due to the proximity of the 
AZO plasma frequency to the NIR (Appendix 0 . On 
the other hand, the Purcell factor for the Ag/InGaAsP 
system is on the order of unity because dk x /du:^0 (Ap¬ 
pendix [A}. Therefore, while Purcell enhancement may 
cause the already stringent gain requirements for the 
AZO/InGaAsP system to increase, they may be consid¬ 
ered a minor concern for the Ag/InGaAsP system. 

To validate the SMM results on MD systems, we per¬ 
formed numerical FDTD simulations (Lumerical®) of 
the transmission through a 10-period Ag/InGaAsP sys¬ 
tem with 30nm layers under the prism coupling con¬ 
figuration and with material parameters identical to 
the SMM. Periodic boundary conditions and perfectly 
matched layers (PML) were enforced in the transverse 
and longitudinal directions, respectively. For each simu¬ 
lation, the center frequency of the finite-bandwidth pulse 
was tuned to match the monochromatic SMM calcu¬ 
lations. Figures HD and [5|b) compare transmission 
through this system at 1500nm and 1550nm, respec¬ 
tively, over a narrow angular range. The SMM is seen 
to consistently agree with FDTD results within a factor 
of two. Over the complete transmission window of the 


(c) 


TEF, lOp 30/30nm Ag/InGaAsP, s p =64, X n =1500nm 




FIG. 5. (a,b) Transmission and (c,d) TEF through prism- 
coupled 10-period 30/30nm Ag/InGaAsP system under in¬ 
version at Ao= (a,c) 1500nm and (b,d) 1550nm. The TMM 
result has been increased lOx in (a,b) for visual comparison. 
(SMM=scattering matrix method, FDTD=fmite-difference 
time-domain, TMM=transfer matrix method) 


Ag/InGaAsP system, agreement between the SMM and 
FDTD worsens (Appendix [B]). This discrepancy between 
the SMM and FDTD may be caused for several reasons, 
including imperfect PML and the finite bandwidth of the 
incident wave used in the FDTD simulations. While the 
exact source of the discrepancy is unknown, we offer sev¬ 
eral possibilities in Appendix [Bj Figure [5ja-b) addition¬ 
ally includes results based on the TMM, which disagrees 
with predictions of SMM and FDTD by several orders of 
magnitude over the entire transmission window (further 
details in Appendix |A| . 


Throughout this analysis we have assumed a steady- 
state gain model, which implies continuous wave external 
pumping conditions. A logical extension of the present 
work is a FDTD study incorporating spatial dependence 
of the gain within each semiconductor layer, as well as 
the use of a self-consistent gain model that takes into 
account transient and nonlinear effects 5 ^ The inhomo- 
geneously broadened lineshape of InGaAsP MQWs how¬ 
ever poses challenges in this regard, as a simple 4-level 
system^ does not capture the electronic density of states 
responsible for the broadening. A potential solution to 
this problem is the use of a superposition of 4-level sys¬ 
tems, each with a slightly different transition frequency. 
The added accuracy of this approach, however, would 
come at the expense of large memory and processing re¬ 
quirements. 
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V. CONCLUSION 

We have analyzed transmission of high-/c waves in the 
presence of loss and gain through metal-InGaAsP sys¬ 
tems. We have shown that, under moderate pumping, 
Ag/InGaAsP may support transmission of high-fc modes 
enhanced by several orders of magnitude, relative to the 
unpumped case, and absolute transmission on the order 
of unity. While the magnitude of the enhancement factor 
for the AZO/InGaAsP system is significant, the absolute 
transmission is rather negligible. We have further shown 
that the Purcell effect does not have a strong deleterious 
role for the Ag/InGaAsP system. The results provide 
justification for experimental efforts on multilayer MD 
systems based on InGaAsP MQWs that circumvent the 
confinement-loss tradeoff fundamental to plasmonic and 


HMMs, as well as a rigorous foundation for further theo¬ 
retical work on the nonlinear dynamics of such systems. 
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Appendix A: Use of Bloch’s Theorem and the TMM 
in analysis of MD systems 

To date, most work on the analysis of MD systems 
and HMMs has relied upon the transfer matrix method 
(TMM). For example, greater-than unity transmission 
through Ag/Ti02 multilayers was predicted via TMM 
by Cortes et a! 5 ^. This result is unphysical because loss 
and gain were omitted. In the presence of loss and gain, 
Savelev et al. used the TMM method to analyze trans¬ 
mission through Ag/PMMA systems^ 7 . While the re¬ 
ported transmission is plausible, no attempt was made 
to validate the results with numerical methods. 

In addition to transmission, the TMM is routinely used 
to calculate the complex amplitude reflectivity of HMMs, 
upon which the photonic density of states and emission 
lifetimes are often calculated. This was done, for exam¬ 
ple, by Galfsky et al. on quantum dots in Ag/A^Os 
multilayerd^, by Lu et al. and Ferrari et al. on dyes in 
Ag/Si system^ 5 ^, and by Cortes et al. for Ag/TiOj^. 
The TMM was also used to elucidate the presence of 
volume plasmon polaritons in AU/AI2O3 multilayers by 


(b) 

log 10 T(k x ,X 0 ), lOp 30/30nm AZO/InGaAsP, Lossless log 10 T(k x A 0 ), lOp 30/30nm Ag/InGaAsP, Lossless 



FIG. 6. Transmission contours for 10-period 30/30nm (a) 
AZO/InGaAsP and (b) Ag/InGaAsP systems with losses and 
gain omitted. The greater-than-unity transmission in the ab¬ 
sence of gain is an unphysical result that motivates the use of 
SMM. 

Zhukovsky et alP^, and to calculate the bandstructure of 
HMMs in the presence of gain 58 . 

While useful qualitatively, the TMM is known 
to become increasingly numerically unstable as the 
structure increases in lengtlF^, as strong loss or 
gain is incorporated^], and as evanescent modes are 
considered^. The essential difference between the TMM 
and SMM lies in the fact that the matrices of the former 
contain exponential functions with arguments of differ¬ 
ent sigEF 8 ^ ] 59 * 6Q [ For large propagation lengths or large 
values of the imaginary component of the wave-vector, 
one of the exponentials will diverge. The matrices of the 
SMM on the other hand contain exponentials with argu¬ 
ments of the same sign, which prevent the instabilities 
that limit the utility of the TMM. Given these facts, it 
is imperative to use a more stable method for calculating 
absolute transmission of high-/c modes in the presence of 
loss and gain, as well as relative transmission enhance¬ 
ment. 

Nevertheless, the TMM does correctly predict the lo¬ 
cation of transmission resonances in k- space. Because it 
is computationally less demanding and simpler to imple¬ 
ment than the SMM, it is useful for studying the disper¬ 
sion of MD systems. Figure [6] shows transmission con¬ 
tours with loss and gain omitted calculated via TMM 
for lOp 30/30nm MD systems. Qualitatively, we observe 
that the AZO/InGaAsP system (Figj6^a)) exhibits little 
dispersion (d\o/dk x ~0), while the Ag/InGaAsP system 
(Fi^6^b)) is highly dispersive (dAo/ dk x ^>0 ). This is a 
direct consequence of the fact that the plasma frequen¬ 
cies of AZO and Ag lie within and outside the NIR, re¬ 
spectively. The greater-than-unity transmission in the 
absence of gain is a quantitatively unphysical result and 
motivates the use of the SMM for absolute transmission 
and relative transmission enhancement calculations. 


Appendix B: Discrepancy between SMM and FDTD 

The discrepancy in absolute transmission between the 
SMM and FDTD calculations is almost negligible at the 
high-& end of the transmission window, but becomes sig¬ 
nificant at small values of k x , in particular at the second 
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lOp 30/30nm Ag/InGaAsP, A 0 =1500nm, N=5e18cm“ 3 



FIG. 7. Transmission through prism-coupled 10- 
period 30/30nm Ag/InGaAsP system under inversion 
at Ao=1500nm. (SMM=scattering matrix method, 
FDTD=hnite-difference time-domain) 


angular resonance of Fig. [7| which is calculated at an 
angular resolution of 7 samples per degree. Potential 
reasons for this discrepancy rest in the finite-size mesh 
of the FDTD simulation and in the finite bandwidth of 
the FDTD source, compared to a monochromatic source 


in the SMM. To check the former potential problem, we 
reduced the mesh size from the default minimum of 4nm 
to O.lnm in both the direction normal and parallel to 
the layer interfaces. We observed negligible change in 
the transmission upon reduction of the mesh size, in¬ 
dicating that a 4nm mesh size was adequate. Given 
that the FDTD is exact once the finite spatial resolu¬ 
tion is taken into account, we suspect that the SMM 
overestimates transmission, particularly at small values 
of k x , due to its monochromatic nature. In the SMM, 
all the energy gained by the input signal must be dis¬ 
tributed across the range of simulated incident angles at 
a single frequency. In a system with optical losses but 
without gain, transmission is generally higher for smaller 
values of k x because the confinement to the lossy metal 
is less at smaller k x . Without being able to distribute 
the energy provided by the gain to other frequencies, the 
added energy in the SMM calculation with gain will in¬ 
crease the angular resonances. For the FDTD however, 
energy gained is distributed across a range of frequen¬ 
cies. In our FDTD simulations, the source has a finite 
bandwidth of 2.67xl0 13 Hz. Arguably, the FDTD result 
is more physical because no source or transmitted sig¬ 
nal is truly monochromatic. Finally it should be empha¬ 
sized that the discrepancy between the SMM and FDTD 
is relatively small compared to the discrepancy between 
the FDTD and TMM. The TMM underestimates trans¬ 
mission by several orders of magnitude compared to the 
SMM and FDTD, and does this consistently over all an¬ 
gles of incidence. 
















